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The PHENIX experiment has recently acquired ~ 1 billion minimum bias Au+Au 
events at yfs = 200AGeV during the year-2004 run. This high statistics data set, coupled 
with a state-of-the-art analysis technique, allows for the extraction of 3D emission sources 
^ . for various particle types. These 3D sources lend fresh insight into the nature of a long- 
range source previously reported by PHENIX. The new results indicate an anisotropic pion 
emission source in the pair center of mass system (PCMS) having an extended space-time 
extent oriented in the outward direction. The two-proton emission source from the same 
data set is essentially isotropic in the PCMS. These results provide a "window" for viewing 
the evolution dynamics of the high energy density nuclear matter created at RHIC. 

V ' 

1. Introduction 

O 

A de-confined phase of nuclear matter is expected to be formed at the high energy 
densities created in Au+Au collisions at RHIC [I]. An observation of the presence (or 
absence) of an emission source of large space-time extent for particle emission can provide 
important constraints for understanding the nature of this phase transition. 

In recent measurements, the PHENIX Collaboration has observed a long range structure 
in the ID two-pion emission source function for Au+Au collisions at -^i = 200 AGeV [2]. 
This long range structure was resolved using the ID Source Imaging technique of Brown 
and Danielewicz |3|4j . The ID technique does not give crucial directional information 
so a detailed analysis of the dynamical origin of this structure was not possible. To 
address this shortcoming, Danielewicz and Pratt introduced the more powerful technique 
of decomposing correlation functions into correlation moments using a cartesian surface- 
spherical harmonics basis |S] • In this representation, each moment (of a particular order) 
corresponds to a specific deformation ( i.e dipole, quadrupole etc ) of the 3D source 
function, and their extraction provide detailed 3D information about the emission source. 

In this paper, we present the first application of the moment decomposition technique 
for the analysis of ti + tt + and pp pairs produced in Au+Au collisions at -Js = 200 AGeV. 

2. Experimental Setup and Data Analysis 

The data presented here were taken by the PHENIX Collaboration during the year- 
2004 run. The colliding beams (y^ = 200 AGeV) were provided by the RHIC accelerator. 
Charged tracks were detected in the two central arms of PHENIX jB] , each of which sub- 
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tends 90 degrees in azimuth and ±0.35 units of pseudo-rapidity. Tracking information was 
provided by a drift chamber followed by two layers of pad chambers. Particle identification 
was performed by an electromagnetic calorimeter and a time-of-flight wall. 

3D correlation functions , C(q), were obtained as the ratio of foreground to background 
distributions in relative momentum q for n + 7v + and pp pairs. Here, half of 

the relative momentum between the two particles in the PCMS frame. The foreground 
distribution was obtained using pairs of particles from the same event and the background 
was obtained by pairing particles from different events. The events used have a z-vertex 
position within ±30cm from the center of the PHENIX spectrometer. Track merging and 
splitting effects were removed by appropriate cuts in the relevant coordinate space on 
both the foreground and background distributions. There was no significant effect on the 
correlation functions due to the momentum resolution of 0.7%. 

In the cartesian harmonic decomposition, the 3D correlation function is expressed as 

C(q) = E E C l a ,.. ai (q) A l ai ^ ai (%) (1) 

l oti...a>i 

where / = 0,1,2,..., «j = x,y or z, A l ai a (fi q ) are cartesian harmonic basis elements 
(fi q is solid angle in q space) and C l a a (q) are cartesian correlation moments given by 

(2) 

The cartesian coordinate system is oriented such that z is parallel to the beam (longitu- 
dinal), x points in the direction of the total momentum of the pair in the PCMS frame 
(outward) and y is perpendicular to the other two axes (sidewards). 

3. Results 

Fig- CD shows the 1=0 moment C° (solid stars) and ID correlation function C(q) (open 
stars) for 7r + 7r + pairs with 0.20 < < 0.36 GeV/c, from Au+Au collisions in the cen- 
trality range 0-30% of total cross section. Here, kx = ( piT + P2T ) is the mean tranverse 
momentum of the two particles in the particle pair. The 1=0 moment is essentially iden- 
tical to the ID correlation function as expected. We conclude that fig. Q serves as a good 
consistency check of the moment calculation procedure. 

Figs. E{a),(b) and (c) show the 1=2 moments C% x , Cy y and C 2 ZZ respectively for tt + ty + 
pairs with 0.20 < k T < 0.36 GeV/c, from Au+Au collisions in the centrality range 0-30%. 
The non-zero values for the 1=2 moments represent source anisotropies of a quadrupole 
nature and point to specific deformations in the direction specified by the moment. More 
specifically, the fact that C% x is negative for all q inv values indicates a negative contribution 
to the correlation function in the x (outward) direction. This is to be compared to the 
results for Cy y and C zz which are both positive and hence represent positive contributions 
to the correlation function in the y (sideward) and z (longitudinal) directions. Thus, the 
overall correlation function, obtained by adding the 1=0 and 1=2 moments, is narrower in 
x and broader in the y and z directions, indicating the 3D source size is larger in x and 
smaller in y and z directions, as compared to the angle-averaged source size. 

Fig- El shows the 1=0 moment C° (solid stars) and ID correlation function C(q) (open 
stars) for pp pairs for 0.4 < k T < 2.0 GeV/c, from Au+Au collisions in the centrality 
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Figure 1. 7r + 7r + 1=0 moment C° and ID cor- 
relation C(qi nv ) for centrality 0-30% and 0.20 < 
kj> < 0.36 GeV/c in Au+Au collisions. 



Figure 2. 7r + 7r + 1=2 moments (a) C xx (b) 
and (c) C 2 ZZ for centrality 0-30% and 0.20 < 
hr < 0.36 GeV/c in Au+Au collisions. 



range 0-90%. Again, the very good agreement between the 1=0 moment and the ID 
correlation function attests to the reliability of the moment calculation procedure. 

Figs.|lja),(b) and (c) show the 1=2 moments C xx , Cy y and C\ z respectively for pp pairs, 
with 0.4 < hr < 2.0 GeV/c, from Au+Au collisions in the centrality range 0-90%. In 
contrast to the n + TT + moments, the 1=2 moments for the pp correlation function are all 
consistent with within statistical fluctuations. Hence, the overall pp correlation function 
is the same in all 3 directions, namely the angle- averaged 1=0 moment C°. 

4. Discussion 

Analysis of the ty + ty + ID correlation function indicates a long range structure in the 
angle-averaged pion source function [2]. The cartesian moments, shown in Fig. |21 indi- 
cate that the correlation function is narrower in the outward direction compared to the 
sideward and longitudinal directions. The consequence of this is a source function which 
is more elongated in the outward direction in the PCMS frame. 

This elongation must necessarily be due to prolonged emissions in the PCMS. A source 
isotropic in its rest frame, but breaking up instantaneously in the PCMS frame, would 
have been shortened in the outward direction by Lorentz factor jboost- On the other hand, 
a source breaking up instantaneously in its own rest frame would lead to a source elon- 
gated in the PCMS by the factor jboost in the outward direction. Here, non-simultaneity 
of emission and movement of the source would lead to overcompensation of the Lorentz 
contraction. The most extreme Lorentz elongation factor would be obtained under the 
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Figure 3. pp 1=0 moment C° and ID 
correlation C(qi nv ) for centrality 0-90% and 
0.4 < kr < 2.0 GeV/c in Au+Au collisions. 



Figure 4. pp 1=2 moments (a) (b) Cy y 
and (c) Cf 2 for centrality 0-90% and 0.4 < 
kx < 2.0 GeV/c in Au+Au collisions. 



assumption of a source freezing out instantaneously in the locally co-moving frame. How- 
ever, such a source would be unlikely. Instead, a more plausible source might be one 
with a radial velocity along the pair momentum, that freezes out over a finite time in its 
frame. Further conclusions await the results of detailed model studies with global and/or 
dynamic parametrization of the emission. 

In contrast to the pions, the cartesian moments for pp pairs, shown in Figs. El and 01 
indicate that the proton source function is isotropic in the PCMS. This suggests a proton 
emission source having a velocity close to that for the pair velocity. 

The current results constitute an important model constraint because a basic require- 
ment will be to explain the pion and proton emission differences and the associated cor- 
relation anisotropies simultaneously. 
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